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Abstract 

We present computing results of the Cherenkov light propagation in air and water from ex- 
tensive air showers developing over the ocean. Limits on zenith angles of the showers, at which 
the registration of flashes of reflected Cherenkov photons by the satellite-based detector TUS is 
possible, are analyzed with consideration for waves on the ocean surface. 

Introduction 

The space experiment TUS is going to be launched in the near future for studying ultrahigh-energy 
cosmic rays (UHECRs) [T|. Its main purpose is to search for UHECRs particles over the whole celestial 
sphere by recording the fluorescence that is generated in developing of extensive air showers (EASs). 
This radiation is isotropic, which allows some photons to reach the satellite-based detector. 

In addition, development of UHE showers is accompanied by the intensive flux of Cherenkov 
radiation, propagating mainly along the shower axis. Nevertheless, when an EAS passes over a high- 
albedo surface the Cherenkov radiation can be detected at high altitudes. An example of this is 
the balloon experiment [2] for recording Cherenkov photons scattered from snow-covered surfaces. 
According to ^ , the total signal of the Cherenkov radiation gives an exact estimate of the primary 
energy, and its spatial distribution enables to measure the shower maximum depth and therefore to 
estimate the mass of the primary particle. 

In the case of the satellite-based detector, recording Cherenkov photons reflected from the ocean 
surface is of great interest as their flash at the end of the fluorescence track would allow one to deter- 
mine the shower trajectory more precisely. Besides, Cherenkov photons scattered in the atmosphere 
produce an additional background. Calculations show that they will be superimposed on the fluores- 
cence signal and blur its maximum. Thus dependence of the "noise" on the shower parameters should 
be found and subtracted from the total signal. Estimates of different components of the Cherenkov 
light, which arise in developing of showers over the ocean, are presented in [4j. 

This paper presents a computing technique for propagating EAS-produced Cherenkov photons 
in air and water. It was used in [3] for estimating the ability of the TUS detector to record the 
reflected component of the Cherenkov light. These calculations may also be useful for experimental 
data processing and for constructing next-generation detectors (see, for instance, the JEM-EUSO 
project [6j). 

1 Calculation technique 

The photon propagation task can be conventionally divided into two parts: 1) generation of the 
photons by electrons of the shower and 2) their subsequent travel in a medium. We apply the Monte- 
Carlo method in both cases, with the destiny of a large number of photons (henceforth, just a photon) 
simultaneously predicted. Each photon is assigned by a certain weight, i.e. the number of photons 
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"survived" at a given calculation stage with respect to the maximum possible number of photons at 
the point of their generation. 

We use the laboratory reference system (LRS) in which the origin O is located at the intersection 
point of the shower axis with the plane Oxy coinciding with the flat ocean surface, and Oz-axis is 
directed to the top of the atmosphere. The LRS and systems to be introduced below are right-handed. 

First we follow the "scheme of photon generation" (section [2]), where the energy and direction of the 
parent electron, the generation depth, direction and weight of the photon are defined. Thus we obtain 
characteristics of the photon at the generation point, which are further used as input parameters in 
simulating its travel. The "scheme of photon propagation in a medium" (sectionjsj) presents a sequence 
of scattering events in which we define the path length of the photon to the next scattering point, 
its direction and weight. Calculations continue until the photon reaches a given level (for instance, 
the height of the satellite) or its weight becomes less than a certain preset value. Reflection and 
refraction of the photons at the "air-water" interface are considered with allowance for waves on the 
ocean surface (section |4]) . 

2 Scheme of photon generation 

When an EAS develops over the ocean surface Cherenkov radiation is generated not only in air. 
Some part of the cascade, depending on depth of the primary interaction as well as energy and direction 
of the shower, penetrates into water and also produces Cherenkov photons. Therefore we consider two 
media: air and water. 

The discussion below will require direction cosines of the shower in the LRS defined in section [T| 
^J'f = - sin ^'sh cos v?sh , fJ-t = - sin O^h sin ip^h , /^f = - cos O^h , 

where the zenith ^sh and azimuth (^sh angles are input parameters. Sought characteristics of a photon, 
namely its coordinates, direction cosines and weight, will be primed. 

We assume that the photon is produced along the shower axis, consequently coordinates of its 
generation point can be calculated from 

,sh ,,sh 



, -hoHX/Xo) (air), , , , > 1^1 
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Here, /io = 8.43 • 10^ cm is the standart atmosphere height, Xq = 1020/ cos g/cm^ is the ocean level 
depth, p„at = lg/cm'^ is the water density. The depth X of the photon generation point is simulated 
from the Gaisser-Hillas formula (see, e.g., [7J) 



N{X) = N, 



max 



X 




^max 


exp 



(Xmax — X)/ A 



with parameters describing the average cascade curve of a shower with primary energy Eq = 10^*^ eV: 
iVmax = 7-1010, Xmax = 800g/cm2, A = 80g/cm2. 

For finding the photon direction in the LRS, it is convenient to introduce local reference systems 
associated with the shower axis (the reference system of the shower, RSS) and with the parent electron 
(the reference system of the electron, RSE). In the RSS the Oz-axis is steered along the shower axis 
and angles 9e and fe characterize the direction of the electron, which, in turn, coincides with the 
Oz-axis of the RSE. The direction of the Cherenkov photon in the RSE is described by angles ^ch and 
ifch- Then its direction cosines /i^, //^ and /i^ in the LRS are given by 

7sh ^^f^^f -7sh 

7sh fJ-f^J-f 7sh Ml 
-l/7sh 
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where 7^= [l — (/i")^] (a = sh, e) and direction cosines of the electron and the photon are defined 
by standart formulae: 

//f = sin0/3 cos(/9^ , = sin 0^ sin , /if = cos 0/3 (/3 = e,ch). 

In calculations we distinguish between "forward-current" (0e^8O°) and "back-current" (0e>8O°) 
electrons, with the latter taken into account only in generating in water. Zenith angles of the back- 
current electrons are distributed isotropically, while that of the forward-current electrons are simulated 
from the appropriate distribution suggested in [8] . Azimuth angles of electrons of both types are 
uniformly distributed. 

For photons, azimuth angles obey the uniform distribution and zenith angles are given by the 
Tamm-Frank formula as a function of the electron energy E: 

sm 



[E/meC^Y - 1 ' 

Here, -Eth = meC^(l — n~'^)~^/'^ is the threshold generation energy of Cherenkov photons by an electron 
in a medium with refractive index n. To determine the energy of the forward-current electrons we 
use the suitable distribution from [8] . The energy of the back-current electrons is found from the 
distribution _ 

N(> E) r E "'"^■■^ 



iV(> ^th) 



0.242 



-E'th -nieC^ 

which was specially derived by means of the Geant4 code [9j (the calculations were performed with 
a lower value of Eq). The number N of the back-current electrons with energies E > Eth is ~ 2.2 % 
of their total number N. 

Finally, the photon is assigned an initial value of its weight. According to the Tamm-Frank theory, 
the number of photons with wavelength A, dq/dX, produced per unit depth, depends on the energy E 
and depth X of the parent electron and is given by 

^^^'^^ = ^^^sin20,h(i^,i?th(X)). (1) 

Therefore the number of the photons in a "packet", the fate of which is simulated, varies from one 
iteration to another. We assume here that the unit-weight photon is generated by an electron with 
infinite energy at the ocean level. Then the weight at any generation point can be defined as 
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We note that a value of P', averaged over a large number of events simulated, allows us to determine 
the total number of Cherenkov photons generated by electrons of an EAS quite easily (see app. [A| . 

Thus, applying this scheme, we are able to specify the coordinates x' , y' , z' , direction cosines fi'^, 
n'y, n'^ and weight P' of a photon at its generation point, which are used further as input parameters 
in simulating its propagation in a medium. 



3 Scheme of photon propagation in a medium 

The brightest lines of fluorescent radiation in the atmosphere lie in the wavelength range A = 
300— 400 nm. So the TUS detector is assumed to be sensitive to photons with these wavelengths and 
we, in turn, restrict our calculations within this range. 

The process of photon propagation in a medium is usually represented as a sequence of scattering 
events, in which the coordinates of a new scattering point and direction of the photon are defined. 
We consider scattering in water and air separately because of differences in their physical properties. 
The input parameters are the coordinates, direction cosines and weight of the photon. 
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3.1 Water 

In simulating the photon travel we take scattering and absorption processes into account. The 
corresponding free path lengths Ig and la are given in [lOj . The density and refractive index of water 
are assumed to be constant: Pwat ~ 1 g/cm"^ and JT^wat ~ 

1.33. 

Let us consider how a single scattering event is simulated. First, we calculate coordinates of the 
scattering point: Xj = + fi'^J {xi = x,y,z). Here, / is the photon path length between succes- 
sive scattering events, which, in the case of a medium with constant density, has the exponential 
distribution: 

p{l)dl = exp{-l/ls)dl/ls. (2) 

The weight of the photon is attenuated due to absorption on the path /: P = P' exp{—l/la)- Then we 
compute new direction cosines of the photon in the LRS: 

7X/^'^ f^'x\ /smOsCosifsX 

I'f^'yf^'z I'f^'x /^y P sin sin for | //'^| < 1, 
f,'J \ cos 9s J 

■ a \ 

sm t/g cos (ps \ 

sin 6*5 sin (/?s for (l — l/x'^l) ^ 10~^, 
cos 9 s / 

where 7'= [l — fi'^~\ The angles 9s and (ps characterize deflection of the photon from its initial 
direction, ips is uniformly distributed and 9s is simulated from the Heneye-Greenstein distribution 
(see, e.g., [TO]): ^ 

p(9s, q) dQs = — -, — ^ ^ ^' — ^ ,^ d9s dips 

'^^ 47r (l + g2 cos 0^)3/2 

with 5 = 0.924. 

Further the next scattering event is simulated. The procedure continues if the condition P^Pmin 
holds or until the photon reaches the ocean surface. 

3.2 Air 

Travelling of a photon in air is considered disregarding its "true" absorption. This assumption is 
quite acceptable in the wavelength range of photons under study here. Besides, we limit to Rayleigh 
scattering by presuming the absence of clouds, aerosols, etc. The atmosphere is supposed to be 
isothermal. 

The coordinates, direction and weight of the photon are determined in the same way as for water. 
However the atmosphere is nonhomogeneous, which prohibits application of ^ for simulating the 
path length to the next scattering point. Nonetheless, within the assumptions indicated, the vertical 
depth |AX„| traversed by the photon has the following simple exponential distribution (see app. [B| ): 

p (AXy) dAXy = exp(-x AX„//i'^) dAX^ . 

Hence we can find the vertical depth and coordinate z of the scattering point: 

X^ = X'^- AX^, z = -ho ln{X^/XM), 

where X^q = 1020 g/cm^ is the vertical depth at the sea level. The photon escapes the atmosphere 
boundary when Xy <0 and reaches the sea level when Xy > 1020 g/cm^. In their turn, the path length 
and other coordinates of the scattering point are calculated as 



I = {z - z')/fi[, X = x' + n'^l, y = y' + fi'yl. 
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The direction cosines of the photon in the LRS are determined from ([3]), with ips distributed 
uniformly and Og simulated from 



p{es) dvts 



IGvr 



(1 + cos^ 9s) sin Os dOg dips ■ 



A new value of the photon weight is calculated with allowance for its attenuation due to scattering 
on the traversed path: 

1 — exp(— X A{,) 

P = P' ^--M-^^^) for/.;<0, 
l-exp[-x(X„o-X;)] 

P = for n'^ ^ 0. 

The weight does not change if the photon leaves the atmosphere boundary or reaches the ocean surface. 
This method of weight calculation makes it possible to "select" photons with low number of scattering 
events. 

Further we simulate the next scattering event. The procedure continues if the condition P^Pmin 
is true or until the photon reaches a required level. 

4 Wavy ocean surface 

We deal with two models of wavy ocean surface: with periodic structure (PS) and with chaotic 
structure (CS). In the first model the surface is represented in the LRS by a motionless wave with 
length Xw and amplitude /i^ along the Ox-axis. The shape of the wave is described by the parametric 
equations of a truncated cycloid: 



t — ho,, sint, 



hw cost, 



where = Xw/2tt and > hw, without dependence on the y-coordinate. In the second model the 
ocean surface is divided into many small flat cells with height rj. The quantities rju = dr]/dx and 
r]c = dr]/dy (the subscripts u and c indicate the upwind and crosswind directions of the wave) follow 
the two-dimensional normal distribution 1 10 



PiVu,Vc) 



1 



exp 



1 



2 / 2 , 2 / ^ 



(4) 



with variances a"}, ■ 



:3.16-10-3ii^, al- 



Uyj depending on the wind velocity Uyj, in m/s. Both 



models assume that photons are reflected and refracted at the interface according to the geometric 
optics laws. Let us consider the calculation procedure by the example of a reflected photon. 

Periodic structure of the ocean surface. Initially we define the projections x' and z' of the 

coordinates x and z = of the point at which the photon crosses the plane Oxy, onto the wavy surface 
with the assumption that one of the wave maxima coincides with the origin of the LRS. Hence the 
x'-coordinate is the residue of x/X^ and z can be derived from the following equations: 



X 

27rr,„ 



r„, arccos(z'//iu;) — y^/i^ — z' ^ for x' ^ Att,/2, 
ryj arccos(z' '/K) + ^Jhl-z'^ for X > Atu/2. 



Further we proceed to the reference system Ox'y' z' in which the normal to the wavy surface at the 
point of the photon incidence serves as the Oz'-axis. Let /Ux, l^y and be initial direction cosines of 
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the photon in the LRS. Then its direction cosines in the Ox'y' z' system are defined as 



cos a cos /3 sin a cos (3 — sin (3 

■ sin a cos a 
cos a sin /3 sin a sin (3 cos /? 

where a = for and a = vr for x'>A^/2, and /3 is the angle between axes Oz' and Oz: 

f3 = arccos 



z' 



Finahy we determine the direction cosines v^, Uy, v!, of the photon reflected and appropriate 
reflection coefficient. Transformation to the LRS is performed from 

cos a cos (3 — sin a cos a sin /3 
sin a cos /3 cos a sin a sin f3 
— sin /3 cos (3 









3) 















Chaotic structure of the ocean surface. At the first step we simulate the tangents 'qu and r/c 
of the angles between a cell and Ox and Oy axes, respectively, from Q and determine the direction 
cosines of the normal in the LRS at each point of the photon incidence: 

Q = ^ ^ ^ ^ 1 

Further the direction cosines of the incident photon in the system of the normal are calculated: 

IJ-'x\ hCxCz iQyCz -iM hx\ 

^J''y\=\-lCy iCx X h/j^ , (5) 

IJ^'zl \ Cx Cy Cz I \f^z/ 

where 7= (1 — Cz)~^^'^- The direction cosines of the reflected photon and its reflection coefficient are 
found from the laws of geometric optics. 

Finally transformation to the LRS is implemented by means of the rotation matrix inverse to the 
matrix from ([5]). 

5 Recording Cherenkov light signal by the TUS detector 

The simulation algorithm presented above was applied in for estimating the signal from reflected 
EAS Cherenkov light by the TUS detector with field of view of Ri9.5°. Within the model of periodic 
structure of a wave, with amplitude /i^ = 0.7m and length A^ = 40m, the registration of flashes of the 
reflected photons was found to be possible for showers with zenith angles < 20° . The wave parameters 
indicated were chosen on the basis of the following data. Satellite observations show that, on 
average, winds under 8m/s in speeds predominate over the ocean. According to the Beaufort scale, 
they are matched to waves with heights up to 2 m. It is known that there exists a wavelength spectrum 
at a given amplitude hw Using a two-dimensional distribution over half-lengths and amplitudes, 
suggested in [12], we chose the mean wavelength, A^ = 0.7m, which corresponds to hw = 0.7 m. 

In this section we compare the calculation results, obtained within both models. Note that the CS 
model seems to be more realistic as constructing a kind of "ripple" on the ocean surface. Moreover, 
it depends on only one parameter, the wind velocity. However the situation where a grid of chaotic 
"flutters" is superimposed on an ordered structure, appeared to us to be more appropriate in the case 
of intense waves. 
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Let us adopt that the wave parameters = 0.7 m and A = 40 m in the PS model are matched 
to the wind velocity = 6m/s in the CS model. Fig. [l] presents angular distributions of reflected 
Cherenkov photons near the ocean surface for vertical showers. Photons with zenith 9 and azimuth 
(f angles impinge on proper cells with area rArAip, where Ar = 0.05 and Aip = it/18. A value of r 
varies from to 1 and specifies sin^. The z-coordinate specifies the decimal logarithm of the photon 
relative surface density. As illustrated, values of 9 in the CS model lies in a much wider range. 




Figure 1: Angular distribution of reflected Cherenkov photons near the ocean surface for showers with angles 
^sh = 0° and ifish — 0° , obtained within the models of a) periodic and b) chaotic structures of the ocean surface. 
The Oz-axis specifies the relative density of the photons on logarithmic scale. 

The restriction on EAS zenith angles, 9sh ^ 20°, has been obtained within the PS model with 
allowance to the condition ipsh = 90°, when the shower axis is in one plane with the wavefront and 
therefore more photons are reflected from a practically flat surface. This is also true for the CS model 
because ac<cru- The azimuthally average number of the photons that will hit the mirror of the TUS 
detector is shown in Fig. [2] for showers with 6'sh = 25°, </?sh = 0° and (/?sh = 90°. 
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Figure 2: Number of reflected Cherenkov photons that hit the detector mirror, as a function of its location R, 
in the CS model for showers with angles 9sh — 25° , (psh — 0° and (psh — 90° . The point i? = at the height of 
400 km is projected onto the center of the LRS. 

Time resolution of the TUS detector Ifis) should be taken into account in order to estimate the 
observed Cherenkov signal value. For small zenith angles of showers, the photons can be considered 
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17.7 ±5.0 


58.5 ±5.1 


92.0 ±5.0 


107.6 ±5.0 


4 


21.5 ±0.6 


62.2 ±0.7 


97.9 ±0.7 


112.6 ±0.7 


8 


21.6 ±0.3 


62.8 ±0.4 


95.3 ±0.4 


113.5 ±0.4 


12 


22.2 ±0.2 


63.3 ±0.2 


96.1 ±0.3 


113.7 ±0.3 


16 


23.0 ±0.2 


65.3 ±0.2 


96.5 ±0.2 


114.0 ±0.2 


20 


24.2 ±0.1 


66.2 ±0.2 


97.1 ±0.2 


114.1 ±0.2 


24 


25.3 ±0.1 


68.0 ±0.1 


99.4 ±0.2 


115.3 ±0.2 


28 


27.2 ±0.1 


69.2 ±0.1 


99.3 ±0.1 


115.1 ±0.1 



Table 1: Number of reflected Cherenkov photons that impinge upon the mirror of the TUS detector in the CS 
model for showers with angles 9sh ~ 25° and (psh = 90° . 



to arrive at the detector mirror simultaneously. However, the signal duration grows starting from 
^sh = 25° — 30° and therefore the signal amplitude reduces even despite of the fact that photons grow 
in number with an increase of the shower zenith angle. Our calculations show that in the CS model, 
with Uu, = 6m/s, flashes of reflected Cherenkov photons at the end of the fluorescence track can be 
recorded by the TUS detector provided that 0shi$25°. 

With a change in wind velocity the restriction on will also change because of proportionality 
of the variances and cj^ to the wind velocity. In tab.jl] are listed the results of the calculations, 
obtained within the CS model at different values of for showers propagating at angles B^h = 25° 
and (f = 90° . The entries of the columns indicate the number of reflected Cherenkov photons that fall 
within the field of view of the TUS detector, as a function of R. As evident, the photons are quite 
small in number at M«, = 2m/s so the maximum value of the shower zenith angle should be decreased. 



Conclusions 

In this paper we have presented the simulation algorithm for generating and propagating EAS 
Cherenkov photons over the ocean. We have considered two different models for describing the wavy 
ocean surface and imposed restrictions on zenith angles of the showers, at which the TUS detector 
would be able to identify the flashes of reflected Cherenkov photons. Our calculations have shown 
that despite the difference in number of the photons impinging on the detector, both models provide 
approximately the same results. 
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A Number of Cherenkov photons in a shower 

The number of photons with wavelength A, produced in a layer {X2—X1), is written as 

X2 00 

i^=/.X'iV(XO/.£'^(X'.£')^(X'.E'). (6) 

The distributions for Cherenkov photons, dq/dX, and for electrons, N{x) and dN{X,E)/dE, are 
introduced in section [21 
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We recall that every electron is assigned not a packet but rather a single photon with weight 

2 



PiX,E) 



Pmax 



sm fch 



sm ( 



ch 



Therefore, ^ can be rewritten as 



d\ A2 



P 



27ra sin^ei^'g'^^ 

Pre, 



X2 



dX' N{X). 



The approximation sign takes in account the fact that ( X^i^i ^i] 1^ tends to P with an increase in 
the number of events simulated. 



B Distribution for vertical depth 

Let us consider a change in intensity of a light beam on a path length dl: 

dl = -Idl/ls. 



(7) 



For Rayleigh scattering, the mean free path length Ig depends on the wavelength A of radiation and 
density p of the medium: 

A4 



L oc 



P 



Introducing the quantity x = {Isp) ^ and using the following relationship between the vertical 
depth Xy of the atmosphere and the height h above the sea level: 



X, 



00 

,{h) = jp{h')dh\ 



we integrate Q over a layer with thickness (/12 — hi): 



In 



I{hi 



h2 

J p{h') dh' 



-X 



X,{hi) - X,{h2] 



-X AX,. 



(8) 



If the beam is deflected at an angle 9 from the vertical direction, ([8]) transforms to 



In 



-xjp{l')dl = —X Xy{hi) — Xy{h2) / cosO = —X /\Xv/ pz , 
h 



where {I2 — h) = (^2 — hi)/ pz- We implicitly assume here the case of the plane-parallel atmosphere, 
which is in reality valid up to ~ 70° . In our calculations photons with /i^ ~ are removed by setting 
their weight to zero. 

Thus the probability that a photon traverses the vertical depth |AX^| without scattering obeys 
the following exponential law: 



p (AX^) d/^Xy = — exp(-x AX^Z/i^) dAX„ . 

Pz 
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